A density functional theory study of electric potential saturation: planar geometry 
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We investigate the possibility of electrostatic potential saturation, which may lead to the phe- 
nomenon of effective charge saturation. The system under study is a uniformly charged infinite 
plane immersed in an arbitrary electrolyte made up of several micro-species. To describe the elec- 
tric double layer, we use a generic density functional theory in which the local micro-ionic density 
profiles are arbitrary functions of the local electrostatic potential. A necessary and sufficient con- 
dition is obtained for saturation, whereby the electrostatic potential created by the plane becomes 
independent of its bare charge, provided the latter is large enough. 
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I. INTRODUCTION 

The importance of effective charge in colloidal or poly- 
electrolyte suspensions has been recognized for some 
time [llllllllllllllli. The electric potential far 
from a colloid immersed in an electrolyte defining the 
inverse screening length k, takes the same form as the 
solution of the linearized Debye-Hiickel equation — say 
ZB^offe~'**-''~^V((l + KR)r) for a charged sphere of ra- 
dius R — but with a prefactor ZeS which is different from 
the bare charge Zbarc of the colloid. In the previous ex- 
pression, r is the distance from the center of the sphere 
of interest, and Ib = q^/{ekBT) is the Bjerrum length, 
defined from the elementary charge q and the permit- 
tivity £ of the solvent (the molecular structure of which 
is neglected). An interesting feature that occurs in the 
framework of (non-linear) Poisson-Boltzmann theory 9] 
is that for highly charged colloids, the effective charge 
Z^fi saturates to a finite value Zgat 0, Interest- 
ingly, this saturation value is an upper bound for effective 
charges found within more refined approaches that incor- 
porate the micro-ionic correlations neglected in Poisson- 
Boltzmann theory |^| . We will come back to the validity 
of the latter approach in the concluding section. 

Since Poisson-Boltzmann is a mean field theory, there 
has been several proposals to go beyond that approxima- 
tion and try to include effects like micro-ions excluded 
volume or electrostatic correlations between micro- 
ions 13j. A natural question arises. Under which condi- 
tions a given theory will account for the saturation phe- 
nomenon? We will not give here a definite answer to this 
question in its full generality. However we will consider 
the special case of local theories — theories which provide 
a local relationship between the density of micro-ions and 



the electric potential — and give a necessary and sufficient 
condition under which the electrostatic potential of an in- 
finite charged plane immersed in an electrolyte will satu- 
rate when the bare surface charge density of the plane di- 
verges. As explained below, the phenomenon of potential 
saturation is slightly more general than that of effective 
charge saturation: effective charge saturation implies po- 
tential saturation, whereas potential saturation does not 
necessarily allows to define an effective charge. 

This paper is organized as follows. We will present in 
Sec. m the framework and the class of theories that will 
be considered. In Sec. IIIII we will formally integrate the 
generalization of Poisson-Boltzmann equation that these 
local theories yield. We will discuss in Sec. IIVI the nec- 
essary conditions to have the potential saturation effect. 
Finally, in Scc.0 we will apply our results to some recent 
theories (1^ that have been proposed to go beyond 
the description of Poisson-Boltzmann theory. 



II. FRAMEWORK 

Let us consider an infinite plane located at a; = with 
a surface charge density a. The half-space a; > is filled 
with an electrolyte made up of several species of micro- 
ions with charges qa and local densities na{x). We use 
Greek indexes for the different species of micro-ions. The 
system is in thermal equilibrium at an inverse tempera- 
ture (3 = {kBT)~^. Without loss of generality we suppose 
cr > 0. 

The average electric potential ili{x) at a distance x from 
the charged plane is related to the densities of micro-ions 
by Poisson equation 



47r 



p{x) 



(2.1) 
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The prime denotes differentiation with respect to x and 
p{x) = Qanaix) + aS{x) is the total charge density. 
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Furthermore the electric potential satisfies the boundary 
conditions 



^'(0+) 

-0(00) - 0, 



= —Ana/e 
^'(00) = 0, 



(2.2a) 

V'"(oo) = (2.2b) 



Suppose now that an approximate local theory is pro- 
vided. This theory gives local relations between the den- 
sities and the electric potential 



(2.3) 



which leads to a closed equation for the electric potential 
when we replace Eq. (|2.3() into Eq. 1)2. 1() 

ii^"{x) = <lc9cc[qa^{x)], for a; > (2.4) 



When ga{u) = exp(— /3m), with the bulk density of 
the species a, we recover Poisson-Boltzmann equation. 
Equation (|2.4I) may then be considered as a generalized 
Poisson-Boltzmann equation. 

The local relation (|2.3I) between the density and the 
electric potential may be obtained in the framework of 
the density functional theory (DFT) by a local density 
approximation (LDA). In this framework, the free energy 
functional is given by 



f({na{x)})dx 



(2.5) 



00 /'OO 



JO 



p{x)Q{x, x')p{x') dx dx' 



where Q is (— 47r/e) times the one-dimensional Laplacian 
Green function with the appropriate (Neumann) bound- 
ary conditions. We have V' (a;) = /g G{x,x')p{x') dx' . 

The minimization of the functional H2.5|) subject to the 
conservation of the total number of particles (controlled 
by Lagrange multipliers fia, the chemical potentials of 
micro-species) gives 



9a V' 



dUa 



(2.6) 



We assume that upon inverting these relations, one ob- 
tains a set of relations of the form H2.3|l between each 
density and the electric potential. Note that Eq. H2.3|l is 
really a set of equations giving, for each a, an explicit re- 
lation between na{x) and the electric potential ipix)- In 
particular each relation involves only one density na{x). 
It can, however, include the charg es and chemical po- 
tentials /i-y of the other particles but not the other 
local densities nj{x). 

Poisson-Boltzmann theory is recovered when the local 
part of the free energy is given by the ideal gas contri- 
bution /({ria}) = /id({n-a}) = Y^^nailniuaA-) - 1], 
where A is an irrelevant length (the de Broglie wave- 
length) . 



The relations (|2.3() should obey a certain number of 
physical constraints, for instance the local charge den- 
sity p{x) should vanish when x 00. We will impose 
the following constraint to the functions ga appearing in 
Eq. ^ 



Yqaga[qa1p{x)]lp'{x) > 



(2.7) 



This condition will be used several times in the sub- 
sequent analysis. In the appendix ^ we show that, 
in the framework of the DFT, this condition is a con- 
sequence of the stability of the system, namely that 
f {{n-a}) / duadn^ is positive definite. 



III. FORMAL INTEGRATION OF THE 
GENERALIZED POISSON-BOLTZMANN 
EQUATION 



Multiplying Eq. by ^'{x), we get 



d r 
dx 



{^'{x)r 



— Yqagaiqai') 



dx 



(3.1) 



which allows for a first integration 



i;'{x) 



— (3-2) 



where we defined 



ga{u) du 



(3.3) 



and we have used the boundary condition l|2.2b(l at a; 
OO. The choice of the minus sign in the r.h.s. of Eq. H3.2|l 
is dictated by the fact that ct > 0. 

Note that condition (|2.7() ensures that the term under 
the square root sign in Eq. (|3.2|l is positive. Indeed, 1)2. 7|) 
implies that d{^' {x)^) / dx < 0, therefore ^i' [x^ is a de- 
creasing function and never vanishes for x finite. Then, 
il^'{x) never changes sign and since cr > 0, ip'{x) < 
and we conclude that ^'{x) is monotonous and increas- 
ing. The electric field —ip'{x) is monotonously decreas- 
ing, always positive and non-vanishing for x finite. This 
remark is in fact a consequence for this particular geom- 
etry of the general proof on the absence of over-charging 
whereby electric double layers are described by a local 
density functional theory of the form H2.5|l [l^ . 

Upon formally integrating Eq. (|3.2I) . we obtain 



F{ij)^-{x + xo) 



(3.4) 



with xq a constant of integration and the indefinite inte- 
gral 



F(^) = 



y/ed(j) 



(3.5) 
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The solution for the electric potential is given by invert- 
ing relation H3.4|l 

4>{x)^F-'[-{x + xo)] (3.6) 

The function F{ip) introduced in Eq. (|3.5|l has a few 
useful properties. It is the integral of a positive quantity, 
so that it is a strictly increasing function of ip for ip S 
[0, oo[^ that can be inverted: exists and it is also an 
increasing function. As a consequence "0(2^) = F~^[—{x+ 
Xq)] is a decreasing function of x. 

The constant of integration xo is related to the sur- 
face charge density a of the plane and is determined by 
the boundary condition (|2.2a(l . An interesting feature 
is that this constant of integration comes as an addi- 
tive offset for the position x. A change in a (therefore 
in the constant of integration xq) results in a transla- 
tion of the curve il){x)-x along the a;-axis. This allows 
a graphical determination of xq: plotting the function 
F~^{—x), the origin a; = of the a;-axis should be such 
that dF~^{—x)/dx ~ —Anaje at this new origin, thus 
satisfying the boundary condition (|2.2a|) . Again, the sta- 
bility condition (|2.7I) ensures that there is a unique solu- 
tion for since (x) is monotonous. 

IV. ELECTROSTATIC POTENTIAL 
SATURATION 

We will say that there is a saturation of the electro- 
static potential if it is possible to have u — s- -l-oo with 
a finite solution for '^{x) for all a; ^ 0^. We emphasize 
that the notion of potential saturation is more general 
than that of effective charge saturation: to define an ef- 
fective charge, one needs to show that the ionic profiles 
behave far from the wall as they would within a linear 
theory such as Debye-Hiickel. The effective charge is 
then defined from the far field created by the charged 
object, as that required within a linear theory to ob- 
tain the same potential at large distances. Our analysis 
does not require such a limitation. However, in the case 
where the local theory is formulated in the framework 
of the DFT we show in appendix 151 that, far from the 
charged wall, the theory reproduces the far field of the 
linear Debye-Hiickel theory, provided that the Hessian 
matrix {d"^ j j duadn^ is positive definite. 



^ We suppose that the function -F(i/') exists for all values of > 0. 
This will be the case if all the functions Ga{qa^) are defined for 
any value of tit. 

^ It is important to exclude x = 0. Indeed, the limit a +oo im- 
plies that V'CO) = — oo but since both ip and are monotonous 
and F{tp) is supposed to exist for all values of i/) > 0, we also 
have VCO) = -|-oo. The case of V'(O) > finite while V'(0) = +00 
is ruled out by the conditions stated in footnote if this was 
the case it would mean that the integral IH.5I defining F(ip) 
would not be defined for values of 1/1 > V(0) which implies that 
GaiqailJa) diverges to —00 for ?/) = i/'(0) and is not defined 
for V > -0(0). 



The function i^('0) is strictly increasing in the interval 
[0, -|-oo[. Therefore there are only two possibilities for 
the behavior of F(ip) when ip — > -l-oo which, as we show 
below, distinguish between the cases of saturation and 
non-saturation. 

If lim F(tp) < +00, there is saturation. (4.1a) 
If lim F{t/j) = +00, there is no saturation. (4.1b) 

'0 — *-|-oo 

Let us consider first the case H4.1a|l . Let 
lim^^oo Fiip) = —Xao < -l-oo. The functions F{ip) and 
il]{x) are sketched in Figure ^ For any finite value of 
cr > the determination of the constant of integration 
Xq gives Xq > x^o and as a increases, xq approaches x^o- 
The case a — -l-oo corresponds to the choice xq — Xoo- 
In the graphical way of determining the constant of in- 
tegration explained before, this means that the origin of 
the X-axis is chosen above Xoo if < cr < -l-oo and it 
approaches Xoo as cr — s- -l-oo. We have clearly the poten- 
tial saturation phenomenon since for cr = -l-oo, where the 
origin of the x-axis is precisely at Xoo, we have a finite 
solution for the potential ij^ix) for any value of a; > 0. 

In the second case l|4.1b|) . Xao recedes to —00. Again, 
for cr > finite, the boundary condition H2.2a|l yields 
a value for xq > Xoo = —00, and as cr increases, a;o 
decreases. We actually have limo-^+cjo xq = —00, the 
constant of integration should decrease to minus infin- 
ity to satisfy the boundary condition. But since ip{x) = 
F^^{—{x + xoj) and i^~^(-|-oo) = -l-oo then the electric 
potential ipix) is infinite for all values of x. There is no 
saturation of the potential. 

In practice, for a general theory formulated in the 
framework of a local DFT, it is not always easy to invert 
Eq. (|2.ti|) in order to have an explicit expression (|2.3|) of 
the density as a function of the potential, and even more 
difficult to compute the function Fiip) to study its limit 
as V' 00. However, there is no need to have explicitly 
the functions ga of Eq. (|2.3|l to discuss the possibility 
of potential or charge saturation. Indeed, as shown be- 
low, we only need to know the behavior of the functions 

gaiqalp) a,S Tp ^ -l-oo. 

First we should treat the particular situation where it 
is possible to have from Eq. (|2.3|) or Eq. (|2.t)|) an infi- 
nite potential ijj +00 while all the densities na are 
finite. If this is the case, then the functions Ga{qa4') 
behave as ip — more precisely they are of order ip which 
we denote 0{^) — when ip — > -l-oo. Therefore the inte- 
grand [— Ga{qa^)\~'^^'^ in Eq. (|3.5|) defining the func- 
tion F{ip) behaves as ip~^/^ — more precisely — 
0{\~Ylia^a{qa'^)\~^^'^^ — as i/) ^ -foo and it is not 
integrable. We have lim^^+00 F{\p>) 00. We then con- 
clude that there is no saturation. To summarize, if a 
local theory predicts the possibility of having all densi- 
ties of the micro- ions finite when t/i — s- 00, this theory will 
consequently not account for the phenomenon of charge 
saturation. 
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FIG. 1: An example of a theory in which the charge satura- 
tion phenomenon is possible. The upper figure represents the 
function Fi^l)) defined in Eq. I|8.5^ and in this example it has 
a finite limit when oo: lim^^+oo F(?/>) — ~Xao < +00. 

The lower figure represents the inverse function F~^{—x). 
The plot of the electric potential ^p{x) as a function of x is 
simply obtained choosing a new origin for the 2;-axis, such 
that dF~^{—x)/dx = —A-Ra/e at this new origin. For a given 
value of a, only the part of the curve at the right of this origin 
(full-line) has a physical meaning. When a increases, this new 
origin gets translated to the left to finally reach the position 
of Xaa when a — +00. 



will not saturate. 

The remaining cases where gyiq-yi/j) diverges faster 
than -0 but slower than ip'^ for any v > 1 can always 
be resolved returning to the study of the behavior of 
[—G'f{qjip)]~^^'^ when ^ —> 00 — the dominant part of 
the integrand of F{ip) in Eq. H3.5|l — and determine if it 
is integrable or not.^ 

In the case where the electric potential saturates, the 
effective charge at saturation can be expressed in terms 
of the function F{tp) defined by Eq. H3.5|l . Let us con- 
sider the saturation regime and suppose that one can de- 
fine an effective charge at saturation, i.e. the far field 
created by the plane has the same form as the one 
predicted by the linear Debye-Hiickel theory (see ap- 
pendix for details). This means that if ?/) ^ we 
have — (47r/£) J2a Qaga{qa'4') ~ '^^'0 so that for x —>■ 00, 



S 3r"t ^ X 



(4.2) 



with ip^'^^ related to the effective charge at saturation by 
the relation ip^^^ = 47r(T^g /(ke). At saturation, in the 
graphical way of determining the constant of integration 
xq, the origin is at Xoo- That is Xao — and 



—X = F{i/j) 
with the particular choice 



(4.3) 



(4.4) 



in the lower limit of integration. To be consistent with 
Eq. H4.2() . we must have F{'ip) ~ lii{ip/->p'''^^) for -ip 
0. Therefore we can extract effective charge at saturation 
from 



^sat ^ ^ ijj^ {V'exp[-KF(V')]} 



(4.5) 



Now, in the opposite case where the theory pre- 
dicts that at least one density Ua — the density of 
counterions — diverges when ip ^ 00, the existence of the 
saturation effect depends on how this density diverges as 
ip ^ 00. If several densities diverge, actually what is 
important is the one that diverges the fastest as ip ^ 00, 
say rij. Remembering that n-y{x) = g^yiq-y'ip), suppose 
that g-f{q-fip) diverges as a power law ip'^ . We prove be- 
low that if > 1, then the effect of charge saturation will 
be accounted for by the theory, otherwise it will not be. 
The argument is simple: suppose g-y{q-y'ip) behaves as -0'^ 
when ip ^ 00, then G.y{qjip) behaves as ip'^^^ . The inte- 
grand in Eq. H3.5(l behaves as as 00 and 
is therefore integrable if and only if {v -I- l)/2 > 1, that 
is > 1. It is also clear from this argument that in the 
cases where the behavior of g-y^q-yip) is not a power law of 
^p but diverges faster than ^p'^ for some v > 1 there will 
be a saturation of the electric potential and in the cases 
where it diverges as ip or slower, the electric potential 



V. SOME EXAMPLES 

In this section we apply our results to some local theo- 
ries proposed in the literature. Wc first illustrate the re- 
sults of the preceding section with some simple examples. 
We consider two cases (as benchmarks, since the analytic 
solution of the problem is known): Poisson-Boltzmann 
theory and its linearized counterpart, Debye-Hiickel the- 
ory. 

For Poisson-Boltzmann theory, ga{qa^) ~ 
Ti'^ex'p{—Pqaip), then the integrand in Eq. I|3.5|) 
behaves as exp(— /3|<7o|V'/2) when -0 — > 00 and it is 



^ For example, there could be cases like [—Gj{q-yip)]~^^^ oc 
(tp lmp)~^ (not integrable) in which there will not be a saturation 
of the potential or cases like [-G^(<j^V)]"^''^ « (?/)(lnV)^)"^ 
which is integrable and there will be a saturation. 
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integrable {qo is the charge of the counterions with 
highest valency, which by the way have a density that 
diverges exponentially faster than the potential ip). 
The function F{'>p) has a finite limit when -0 —> oo; we 
therefore recover there the well known fact that effective 
charge saturates when cr — > oo. 

For Debye-Hiickel theory, gaila^) = n^ai^ — Pla"^) 
behaves as tp when -0 — > oo and the integrand in Eq. (|3.5|) 
behaves as ip~^-< it is not integrable: there is no charge 
saturation. 

Now let us turn our attention to some more interest- 



ing examples. Barbosa et al. have proposed a local 
theory to account for counterions correlations in a one- 
component plasma (OCP) description of the electrolyte 
— that is a system of charged counterions immersed in 
an uniform oppositely charged background. This theory, 
referred to as the Debye-Hiickel-Hole-Cavity approach, 
is stable — it satisfies condition H2.7|) — and the density 
n of the counterions is described in the framework H2.5|l 
of the DFT. The local free energy density is of the form 
/("-) = fidin) + /dhhc("-) where /id is the ideal gas part 
of the free energy and 



^/dhhc('t-) _ (ko)^ 



duj 



2(^3 - f ) 

I — 



(i + i7(<z.)i/3)((:^2 + ^ + i) 



(5.1) 



In this equation. 



f^(cD) ^ (cS-l)3 + i-^(w3-l), w = (1+3/bk)1/3 (5.2) 
31bk 



and K = y/A-Klsn is the inverse Debye length. The theory 
has a parameter a which may be interpreted as a sort of 
radius of the micro-ions (non-strictly speaking, because 
the counterion density can actually be higher that a~^). 

The expression (|5.1f) is quite complicated and there is 
no hope to be able to obtain an analytical solution for the 
function F{'ip) or even invert the relationship between the 
potential and the density obtained from the stationary 
equation (|2.6|) . However we only need to investigate the 
limit ij) —> oo. 

First it is straightforward to see that if the density n 
is finite, the potential is finite. Now, if n — > oo we have 
/3/dhhc("') ~ n^/'^ a?{TT /d))'^/'^ and it is the leading term 
in f(ri). The minimization equation H2.6() then yields 
(i\q\'ip ~ n2/3a2(5/3)(7r/6)2/3 for n ^ oo, so that for 
ip ^ oo, we have 



3/2 



3/2 



(5.3) 



The important fact is that the density of counterions be- 
have as ip'^ with = 3/2 > 1. This theory therefore leads 
to a saturation of the electrostatic potential, and more 
precisely to a saturation of effective charge''. 

On the other hand there are theories that do not ac- 
count for the saturation effect. As an example, let us 
consider the approach proposed by Borukhov et al. [T^ . 



* To allow for the definition of an effective charge, Poisson- 
Boltzmann generalized equation 12. 4i should behave as 

0. This is discussed in 



Helmholtz equation (i/i" cx: rp) when 4> ■ 
appendix El 



This is a local theory which incorporates approximately 
steric effects due to volume exclusion between the micro- 
ions of the electrolyte. The local part of the free energy 
is f{{na}) = /id({"-Q}) + fcKc{{na}) with an excess part 
involving the size a of micro-ions 



(5.4) 

This theory is also stable, namely {d^f/dnaduj) is pos- 
itive definite it therefore satisfies the condi- 
tion H2.7|l . From Eq. (|5.4() . there is a higher bound for 



the value of the micro-ions densities which is 



when 



the counterions density approaches a^^ the electric po- 
tential V' diverges. Therefore, according to the discussion 
of last section, this theory does not allow the electrostatic 
potential to saturate. 

Our approach also predicts, in the case of non- 
saturation, if the effective charge will grow faster or 
slower than the bare charge. Let us illustrate this point 
considering again the framework presented in Ref. 
To simplify the argument, we consider a two-component 
electrolyte with counterions/coions of charge —q/q and 
bulk density n^. Poisson's equation then takes the form 



pqiP"{x) = 



sinh(/3(7'0) 



1 - C + C cosh{Pqip) 



(5.5) 



where = 2ni,a^ and k 

In this case where it is possible to have finite counterion 
density when '0 ^ oo and as discussed in the last section, 
the function F{'ip) defined in Eq. H3.5|l behaves as 0^/^ 
as ip diverges. More precisely, we have here 

1/2 

V//^ + C'(l) (5.6) 



■ip—^o 



2^\q\ 



so that 



(5.7) 
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when X — > — oo. On the other hand, one recovers Debye- 
Hiickel theory for smaU potentials. This means that 
F^^{~x) ~ Aexp(— Kx) for x — s- oo with A an arbitrary 
constant since F is determined up to an additive con- 
stant. Suppose that initiaUy the plane has a small bare 
charge density a^^\ For Debye-Hiickel theory the elec- 
tric potential reads tpDHix) = (47rcr(°)/eK) exp(— kz), so 
that eKA/(47r) may be considered as the effective charge 
describing the far field of the plane. As explained before, 
ipix) — F~^[—{x + Xq)] where the constant of integra- 
tion xq can be determined graphically. Let us choose 
A = 47rcr(*'V(^^)) so that plots of F^^{—x) and ipBuix) 
superimpose when a; ^ oo (see Figure (SJ. Now, consider 
a large bare charge density a of the plane. To obtain from 
the curve of F~^{—x) the correct plot of ip{x) we should 
change the origin of the a;-axis such that the value of 
dF~^{—x) I dx = —Aira/e at this new origin. Let us sup- 
pose a large enough so that this new origin is in the region 
x — > — oo where F~^{—x) behaves as x^ . As cr increases, 
we have to shift the origin further to the left or equiva- 
lently to translate the curve of F~^{—x) to the right such 
that at the origin, both curves of F^^{~x) and V'dh have 
the same slope. But since — V'dh(^) behaves as exp(— kx) 
which has a faster increase as x ^ — oo than -F-i {-x) 
— which behaves as —x — it is clear that as a increases 
we need to shift more to right the curve of F~^{—x) than 
the one for -^DHla^)- On the other hand for x — > oo, we 
have ^{x) = Bexp(—Kx) with B oc exp(— kxq)- Since 
the translation of the curve of F~^(—x) becomes larger 
and larger than the one for '0dh(3:), this means that the 
factor B increases faster than the bare charge a. The 
corresponding effective charge therefore increases faster 
than CT as cr is raised. We have also checked this feature 
from a direct numerical solution of Eq. H5.5|l . The above 
argument may be rationalized, as shown in appendix lO 

In the cases of non-saturation, generalizing the argu- 
ment presented before, if the faster divergent counterion 
density, rij, behaves slower than ip when -0 ^ oo, then 
the theory will predict an effective charge that diverges 
faster than the bare charge. If, on the other hand, di- 
verges faster than ip — but slower than ip'^ for any i/ > 1, 
so there is no saturation of the effective charge — then the 
effective charge will diverge slower than the bare charge.^ 

To sum up, a theory such as that of Ref. [l^ cures 
a deficiency of Poisson-Boltzmann theory by setting an 
upper limit for the density of micro-ions, at the expense 
of producing unphysical effective charges when the non- 
linear character of the equation prevails (only for small 
electrostatic couplings would this theory give reasonable 
effective charges, that would then coincide with the bare 
charge of the plane since the theory would reduce to 
Debye-Hiickel). The same would happen for any DFT 
with the hard core part included through a local free cn- 



^ The first example of footnote l3l is a case where the divergence of 
the effective charge will be slower than the bare charge. 
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FIG. 2: Illustration of a theory (Ref. O) where the effec- 
tive charge diverges faster than the bare charge. The dashed 
curve represents the potential 7/'DH(a;) predicted by Debye- 
Hiickel theory and the full line represents the potential tp{x) 
put forward in Ref. [T^ . The potential is measured in units 
of kT/q and the distance x in units of the Debye length 
We have chosen here a'^nt — 1/4, where nj, is the bulk den- 
sity of ions. Figure represents a case where the bare 
charge a is small (actually in the figure 4naf3q/{eK) = 1). 
Both curves have the same behavior as a; oo, namely, 
7/)(a;) = exp(— a;). Figure |2l3 represents a case where a is 
large {AixaPq/ (en) = 10). To obtain Figure |5j3 from Fig- 
ure|2JV, it is necessary to shift the curves by an amount indi- 
cated in Figure |5J\ by the horizontal arrows. Since when 
X ^ —oo the Debye-Hiickel electric field increases expo- 
nentially —il''u-[i{^) ~ exp(— a;) whereas the electric field of 
Ref. jl2l increases only linearly —ip'{x) ~ —x, the curve for 
that last theory should be shifted to the right much more 
that the one for Debye-Hiickel theory. This makes that the 
behavior at a; oo is ip-Dnix) ~ 47rcr/(eK) exp(— a;) and 
7/)(a;) ~ 47r(7cfr/(eK) exp(— x) with an effective charge larger 
than the bare charge GcS > cr- 

ergy density (one may wish for instance to use Carnahan- 
Starling like expressions to account for excluded volume 
p^): such approaches indeed lead to a divergence of the 
free energy density in the vicinity of close packing, so that 
F{'ip) oc V"^^^ for large V'- Since we have to expect that 
effective charges are bounded from above when the bare 
charge becomes large the above analysis seems to 
show that when improving upon Poisson-Boltzmann the- 
ory, incorporation of steric effects through a local func- 
tional and/or neglect of micro- ionic electrostatic correla- 
tions is insufficient. 



VI. SUMMARY AND DISCUSSION 

We have found a necessary and sufficient condition that 
a generic local density functional theory should obey to 
describe the effect of electrostatic potential saturation for 
a charged planar interface. Namely, if in the local rela- 
tion between the potential and the densities it is possible 
to have i/i — > oo and all densities finite, there will be no 
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saturation. On the other hand, if as — > oo the highest 
valency counterions density diverges as Tp'^ for some v > 1 
or faster, the theory will predict a potential saturation. 
In the other case where the increase of the density is as 
or slower, there will not be an electric potential satura- 
tion, hence no charge saturation. For the limiting cases 
where the density diverges faster than ij: but slower than 
for any v > 1 one should consider the general condi- 
tion: [- with Ga defined in Eq. 
should be integrable as V' ^ oo to ensure that the electric 
potential saturates. 

A natural continuation to this work would be to ad- 
dress the same question for a colloidal object of any 
shape. We expect that, for large colloids, the results 
found here will hold. However for small charged ob- 
jects immersed in an electrolyte, the question remains 
open 21]. It would also be interesting to explore in more 
detail the behaviour of effective charges in the framework 
of other theories, for instance non-local theories. 

At this point, it seems appropriate to come back to the 
validity of one of the simplest approaches encompassed 
by our formalism: Poisson-Boltzmann (PB) theory. In 
particular, it seems that the validity of the saturation 
regime within PB theory is sometimes confuse in the lit- 
erature, which we might trace back to Onsager who wrote 
in a seminal paper [ ITiJ 



As soon as the higher order /non- linear/ 
terms in the Poisson-Boltzmann equation be- 
come important, we can no longer expect the 
ionic atmosphere to be additive, and then the 
Poisson-Boltzmann equation itself becomes 
unreliable. 



This statement (irrelevance of PB as soon as non-linear 
behaviour sets in) is justified for simple electrolytes, the 
system studied by Onsager in 0|. However, it turns 
incorrect when there is a large size and charge asym- 
metry between the constituents of a charged mixture, 
which is the case for colloidal suspensions 18]. This 
appears in the analytical work of Netz in the no salt 
limit Indeed, PB neglects correlations between mi- 
croscopic species, that become prevalent at high electro- 
static coup lings, i.e. when F — {q'^al^Y^^ becomes large 
H IIIIM Il2fl]. On the other hand, PB enters the sat- 
uration regime when the coupling between the charged 
plane and the small ions becomes large, which may still 
be compatible with the neglect of micro-ionic correla- 
tions, provided kIb <C 1, as may be seen replacing the 
saturation value dsat oc k/Ib in the definition of F. A sim- 
ilar conclusion is reached enforcing that the correlational 
contribution to the pressure in the bulk of the electrolyte 
(the celebrated term in — k'^/crT' appearing within Debye- 
Hiickel theory of electrolytes |H,f2 2]) should be negligible 
compared to the ideal gas contribution n^kBT: <^ rif, 
imposes kIb ^ 1. 



APPENDIX A 



We prove here that condition (|2.7|) is a consequence of 
the stability of the theory formulated in the framework 
of the local DFT lE^Jl. First, note that 



d 
dx 



^qa9a{qa1p{x))ljj'{a 

a 

= ^ qan'^ip' + ^ qanai^'' 

a a 



a, 7 



— p 

s 



(Ala) 
(Alb) 

(Ale) 



where we have used the minimization equation (|2.6|l and 
Poisson equation (|2.1(l . If [d"^ f / duadn^) is positive defi- 
nite, this latter expression is negative and therefore 



\qc,ga{qafp{x)W{x) 



(A2) 



is a decreasing function of x. Furthermore, the bound- 
ary conditions at a; 00 impose that (|A2I) vanishes at 
infinity, so that Y.a^a9ai.<lai^{x))il}' {x) > 0. 

The condition that f / duaduj is positive definite is 
a sufficient condition to ensure that H2.7() is satisfied. Let 
us remark that this same condition is sufficient for the 
non-existence of like-charge attraction between colloids 
and the absence of over-charging T^. It is also a suffi- 
cient condition for the compressibility of the system to 
be always positive within the cell model [l^. Further- 
more, as will be shown in appendix ^ it is a sufficient 
and necessary condition to recover Debye-Hiickel linear 
theory where the electric potential is small. 



APPENDIX B 

In this appendix we prove that, far from the charged 
wall, any local theory reproduces the results of the linear 
Debye-Hiickel theory provided that the Hessian matrix 
d-^f/duadnj is positive definite. 

Far from the wall, the densities converge to their bulk 
values and the potential vanishes. Then the minimiza- 
tion equation H2.6|) for ■(/; = becomes fia = {df/dna)b, 
where the subscript b means that the partial derivative of 
/ is evaluated at the bulk values of the densities Ua = n^- 
This provides the relationship between the chemical po- 
tentials and the bulk densities. 

For small values of the potential, the difference — 71^ 
is small, then we can expand 



dl_ 

dUa 



dl_ 
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duadn^ 



(Bl) 



and replacing into the minimization equation (|2.t)|) this 
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gives 



duadn. 



-gaV' (B2) 



7/6 



Provided that the Hessian matrix ( a ^ / — | is invertible 
we have 

where ( „ ^ / — | denotes the (a, 7) matrix element of 
the inverse of the Hessian matrix ( I — I • 

y onaOriy j ^ 

Replacing into the generalized Poisson-Boltzmann 
equation H2.4|l we obtain a Debye-Hiickel like equation 



= n'^ipix) (B4) 
with an inverse Debye length defined by 



2 47r^ f d'f 



a, 7 



duadn. 



(B5) 



7/6 



To recover the results from Debye-Hiickel theory, what- 
ever the values of the charge qa might be, it is necessary 
and sufficient that the Hessian matrix ( „ ^ / — ) is pos- 

y onaOn^ J ^ ^ 

itive definite, which ensures the existence of its inverse 
and that > 0. 



APPENDIX C 



that is FbaoH') behaves as -Fdh(V') when i/j ^ up 
to terms that vanish when ?/; — > 0. On the other hand 
we have the behavior H5.6(l for -Fbao when -0 — > cx). 
Let us define 'i/'baro = AiTa/^Ke) and suppose that the 
bare charge a of the plane is large enough such that 
Pqiphnrc ^ 1- For this value cr of the charge density 
of the plane, we have to find the constants of integration 
a;^^° and x}^^ of Eq. ^ for the theory of Ref. [H 
and for Debye-Hiickel theory respectively. We can pro- 
ceed as follows: we first determine the value of ■0DH,bao 



such that i^DH,BA0(V'°"'^^°) 



DH.BAO 
Co 



, then we de- 



DH BAO 

duce Xq ' . For Debye-Hiickel theory, applying the 
boundary condition (|2.2a|l we have 



dip 
dx 



fdFi 



DH 



1p — 1p^ 



(C4) 

Then for = ^baro we have Fdh(^/'°") = -2:^0" and 
therefore — kxq ^ = In [ipharc / j'^^ ^) ■ F or the theory of 
Ref. the boundary condition (|2.2al) yields 



dx 



x=0 



/cLFbaoV' 
2k 



= '^V'baro 



BAO 



(C5) 



^BAO _ 



Therefore we have ip^^'^ ^ C/^'zV'baro/^ and 
KFBAoi^^^°) ~ C/39V'baro. Clearly we have x^^^° > 



,DH 



as Ipha 



00. Furthermore these constants of 



integration manifest in the large- a; behavior of ipix) as 

ij{x) ~ ^(0)g-fta;og-Ka; giving, for X OO, 



In this appendix we prove that the theory presented in 
Ref. predicts an effective charge that diverges expo- 
nentially faster than the bare charge when cr —> 00. 

The function defined by Eq. H3.5|l for the theory of 
Ref. [13 reads 

^^^1/2 rip ^0 

ln[l-C + Ccosh(/3#)] 
while the one for Debye-Hiickel theory is 

FdhW = ^ln(^/0(°)) (C2) 

where in the last equation we have made the particular 
choice -0 = V'*'"'' of the lower bound of integration in the 
definition (|3.5|l of Fdh- Suppose that Pqip'^^^ is chosen of 
order one — say, for example, (iq-ip^^'^ = 1. Now we choose 
the lower bound of integration in Eq. ljCl(l for Fbao such 
that 

FbaoW - K-^nf^/V^")) +0(1) (C3) 



i'mix) = iZ-baroe (C6) 

for Debye-Hiickel theory (as expected) , while 

^BAo(2:) - ^(0)eC/39^b„c+o(Vba,c)g-«x 

X — *oo 

for the theory of Ref. This theory predicts an ef- 

fective charge that diverges exponentially fast when the 
bare charge cr — > cx). 
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